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Abstract 
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1 Introduction 



It is clearly interesting to find the widest possible class of boundary conditions compatible 
with integrability asocciated to a given model. 

Not any boundary condition (b.c.) obeys this requirement. Periodic and twisted (under 
a symmetry of the model) b.c. are usually compatible with the Yang-Baxter equations 
|1|, 0]. In addition, there are the b.c. defined by refiection matrices [1^, p^ . 

These matrices can be interpreted as defining the scattering by the boundaries. In a 
recent publication ||TI] the interpretation of this matrices as boundary S-matrices in two 



dimensional integrable quantum field theories was developed. They also imply boundary 
terms for the spin hamiltonians which can be interpreted as the coupling with magnetic 
fields in the edges of the chain. 

In addition, quantum group invariance arises for specific choices of fixed b.c. (See for 
example 0, |^, for the trigonometric case and |§] for the elhptic one). A quantum 
group-like structure is still is to be found for which Baxter's 8-vertex elliptic matrix 
could act as an intertwiner (for a recent attempt see [Q) giving an affine quantum 
invariance to the infinite spin chain and the boundary terms for the quantum group 
invariance of the finite chain. This program has been done in the elliptic case for the free 
fermionic model, see 0, |^. 

A general setting to find boundary terms compatible with integrability was proposed by 
Sklyanin P] . To find these boundary conditions one has to solve the so called refiection 
equations: 

R{e - ff) [K- (9) ® l]R{e + 6') [K~ {6') ® 1] = 

[K-{e') ® l]Ri9 + 9')[K-{9) ® 1]R{9 - 9') , ^ ' 

R{9 - 9')[1 ® K+{9)]Ri9 + 9')[1 ® K+{9')] = 

[l^K+{9')]R{9 + 9')[l^K-^{9)]R{9-9') , ^ ' 

where R{9) is the R- matrix of the chain and K'^{9) give the boundary terms (see 
below) . 

As is known, the XYZ model is obtained from the elliptic eight-vertex solution of the 
Yang-Baxter equation: 

[1 ® Ri9 - 9')] [R{9) ® 1] [1 ® R{9')] = 
[R{9') ® 1] [1 ® R{9)] [R{9 -9')^1] . 
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The XXZ and XXX models follow respectively from the trigonometric and rational 
limits of this R-matrix. 

We present in this paper the general solutions K^{6) to these equations for the XYZ, 
XXZ and XXX models. We find for the elliptic case two families of solutions, each family 
depending on one continuous and one discrete parameter, see equations (^) and ([361). 
For the trigonometric and rational limit we find a family of solutions depending on four 
continuous parameters, see equations (|4^) and (^2]) respectively. 

We remark that the trigonometric limit of the elliptic solutions of (|T]),(0) does not provide 
all solutions to the trigonometric/ hyperbolic case. 

From these K^{9) solutions we derive the boundary terms in the XYZ hamiltonian wich 
are compatible with integrability. Finally we analyze the relation of the present eight 
vertex results with the general K- matrices of the six- vertex reported in ref . |^ and consider 
in addition the rational limit. 

2 General solution to the reflection equations for the 
Eight Vertex model (XYZ chain) 

The R-matrix for the XYZ chain can be written as ill: 



R{e) 



1 /c sn7 sn6' ^ 

^ sn(6l+7) sn(6l+7) ^ 
n sn6> sn7 r\ 

^ sn(6l+7) sn(e+7) ^ 



(4) 



^ k sn7 sn0 1 y 

where sn (and cn, dn in the formulas below) stand for Jacobi elliptic functions of 
modulus < < 1. 

This solution of the Yang-Baxter equations enjoys the following properties 

a) Regularity: i?(0) = 1. 

b) Parity invariance: PR{9)P = R{e) where P^^ = 6^6^^. 

c) Time reversal invariance: R^^ = R^f. 

d) Crossing unitarity: R{e)R{-e - 2r]) = p{e) 1. 
Where R^^ = R^^, // = 7 and: 

HO) = 1- 



sn2(7 + e) 
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sn(6' + 27) sn6'.^ ,22 1, n\^ 
= ln^(^^ + e) [l-^'^Asn^(7 + g)] . (5) 

From @ and a) unitarity follows : 

R{e)R{-e) = p{e)i (6) 

sn^7 — 

sn(7 + 6) sn(7 - 6) ' ^ ' 

It is shown in |^ that when the R- matrix enjoys properties b),c),d) and (P) we can 
look for solutions to equations (|1|) and in order to find open boundary conditions 
compatible with integrability . 

Since b) holds, equations (0) and (||) are equivalent. We now look for the general solution 
of these equations in the form: 

Inserting equations (||) and (||) in (|l]) we find twelve independent equations: 

h^yz' + c'^d^ zz' = c^d^yy' + h^zy' (9) 
d~vv' + d'^xv' = d'^vx' + d^xx' (10) 
h~yz' + c~d'^zz' = c~d^yy' + h^zy' (11) 
b^c^f f ' + b~c^xf ' = c^b~fa;' + c"b'''xa;' (12) 
c^yx' + h'^d^zx' + d~vz' + d~^xz' = 

c''yx' + h~^c~vy' + h^c^xy' + h^d'^zx' (13) 
h^yv' + d^d'^vy' + xy' + c~^d^zv' = 

h^yx' + h^h^vy' + c^c^xy' + c~d^2;x' (14) 
h^d'^yx' + qT zx' + b~*'c~f 2;' + b~c"^x2;' = 

h^d^yx' + d~vy' + d~^xy' + c'''za;' (15) 
h~yv' + c~c^vy' + b~b^x?/' + c~d'''2;f ' = 

b"''?/x' + vy' + d~d''"x?/' + c"^d~za;' (16) 
c'd'^ya;' + h~ zx' + b^b^f 2;' + c~c^xz' = 
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c^d yv' + h'^zv' + z'x + d^d vz' 
c~yv' + h~c^vy' + h^c^xy' + b^d^ 



zv = 



(17) 



c^yv' + b''"d 2;f ' + d^f 2;' + d xz' 



(18) 



c'''d yx' + b'^zx' + f 2;' + d d^xz' = 



c d'^yv' + b 2;f ' + c c^vz + b b"'"xz' 



(19) 



b^d yv' + d^tiy' + d xy' + c'^^t'' = 



b d'^yv' + c zv' + c^b f + b^c xz' 



(20) 



where: 



i?(e ± e') 



/ 1 d± \ 

b± c± 

c± b± 

V d± 1 / 



(21) 



and x' = x{0'), y' = y{0'), etc. 
We start by assuming one of the elements of K in equation (P) is equal to zero. There 
will be four cases depending on which element is zero, but only two of them turn out to 
be different: 
a)x = = x' 

Using equation (0) we have v = = v' and we are left just with equations ([Tl|) and (||) 
as independent equations. In order that these two equations be satisfied it must be that: 



where from now on an arbitrary multiplicative function of 6 will be omitted. 
The case where v = = v' is equivalent to this. 

h)z = = z' 

From eq.(^ y = = y' and from eqs. (P^) and (PUf): 




(22) 



b+ + c+d z/y 



which implies {z/y^ = 1- Two solutions are then obtained: 




(23) 
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c"'"b v/x + c b"*" 

(24) 



b+c t'/x + b^c+ 



d^v/x + d+ 



which imphes [v/xY = 1, and then: 




K{0) = K , . • (25) 



The case where y = = ?/' is equivalent to this. 
We now assume x{9) ^ and y{9) ^ 0. Then equations ([TT|) and (^ imply that: 

^(^) = ±1/(0) 7^ , (26) 

and (|12|) and ([1^) require: 

f (0) = ±x(0) ^ . (27) 
The matrices K{9) in this case have the form: 

y e22;(0) e3w(0) y 

where e\^e\,e\ = 1. Omitting an arbitrary multiplicative function of 6 we have only 
eight different possibilities: 

a) ei = ea = 63 = 1 

b) ei = €2 = 1 and 63 = —1 

c) ei = e3 = 1 and €2 = —1 

d) 62 = £3 = 1 and ei = —1 

e) ei = 1 and 62 = £3 = —1 

f) €2 = 1 and ei = €3 = —1 

g) 63 = 1 and ei = €2 = -1 

h) ei = €2 = €3 = -1 

Inserting ( [28| ) in the rest of the equations, one finds only two different equations for 
w{9) = y{6)/x{6) in all cases. They are: 
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(a\i (a'\ esb+c + b c+ - eie2d+ - eiesegd 

wiO) wiO) = r— I r— , (29) 

^ ^' ^ ^ C+ + eicab+d- - - eie2b-d+ ^ ^ 



, b b+ + e^c c+ — d d+ — 63 , , 

wie) w{9') = — --f . 30 

^ ^ ^ b+ + eie2C+d- - egb" - eie2e3C-d+ ^ ^ 

For the previous equations to have a solution, the r.h.s. of (^9]) and (|30D must be 
identical. This can be seen, with some work, to occur for all cases a)-h). One can see also 
that these expressions factorize as: 

w{e)/w{e')= '""l I ^"f , (31) 

for the cases where €3 = 1, and as: 

w{e) w[e) = -— — 5tt/t— r — ^ , 32 

for cases where = —1. We therefore have a 6'- independent free parameter in the 
general solution that we call A. The solution then reads: 

, sn0 , , 

= ^T^ ^ , 33 

^ ^ [1 + eie2k sn^e] ^ ' 

when £3 = 1, and: 

- cn6' dn6l 

w(Q) = At- — - , (34) 

^ ^ [1 + eie2k snW] ^ ' 

when 63 = —1. 

It can be noticed that these solutions are easily obtained by the residue of (^) when 

^ if 63 = 1, or the limit 6' ^ of the same equation when €3 = —1. 
We summarize the general solution of the factorization equations for the 8-vertex model 
as: 



Xha snO [1 + ek sn^O] 



ka{9) = { .;:^"r L i , (35) 



and: 
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[1 + eksn^e] eA^^cn^dn^ 
' AHBcn^dn^ -[1 + ek sn^O] ' ^ ^ 

where = 1 and Aj^/a, A/^b are arbitrary parameters. That is, we find two famihes of 
solutions each one depending on a continuous and on a discrete parameter. (The discrete 
parameter takes only two values). 

These solutions lead in the trigonometric limit A; = to only some specific cases of the 

general solution for the six vertex R-matrix as discussed in the next section. 

We now look for the hamiltonians obtained by the first derivative of the transfer matrix 



H = E + i(AT(0)-')AT(0) + . (37) 



n=l 

where: 



tr[K+{-ri)] 



(0) , (38) 

and the term [K^ {Q)~^)Ki {Q) generalizes the formula for the hamiltonian of Sklyanin 
to the case when -ft'^(O) ^ 1. This formula is only defined when tr[K^{—ri)\ ^ and 
det[K-{Q)] ^ 0. 

We see in equation ( p6D that for the second family of solutions the trace of K is zero. 
For this second family we will then not have a well defined hamiltonian from the first 
derivative of the transfer matrix. 
When 



tr[K^{-r])] = , (39) 

and: 

tro[K+\-ri)hm] (X I , (40) 

a well defined hamiltonian with only nearest neighbours interactions is obtained from 
the second derivative of the transfer matrix as shown in But for the present second 
familiy of solutions the condition (^) does not hold. Furthermore -^'^(O) = wich gives 
only a trivial boundary term at the left end. The same happens with solutions (|23| ) , (|25|) 
where one of the elements is zero. 
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If the condition (^) holds but not eg . (|40|) , one obtains from the second derivative of the 
transfer matrix a hamiltonian with terms that couple every pair of sites in the bulk with 
the boundary. That is, a non local hamiltonian arises. 

The hamiltonians associated to the first family of solutions (^) are given by: 



N-l 



(41) 



i=l 



Here a and (3 can take the values x or ?/ in all possible combinations and the ^± are 
arbitrary parameters proportional to Xha- 

As is clear, by rotating the axis, we can make the indices a and (3 in equation (^) take 
also the value z. 

Equation ( PD gives the most general choice of boundary conditions compatible with 
integrability for the XYZ chain besides periodic and twisted boundary conditions. By 
twisted boundary conditions we mean: 

a%^^ = Ma^M-^ . (42) 

Where a = x,y, z and the twisting matrix M stands for a discrete symmetry of the 
eight-vertex model. That is, M = cr^ or cr^. 

In conclusion, the XYZ hamiltonian is integrable with boundary conditions that corre- 
spond to the coupling with a magnetic field on the end sites oriented along parallel or 
orthogonal directions. 



3 General K-matrices for the Six Vertex model (XXZ 
and XXX chains) 

In this section, we briefly review the results of concerning the general solution for the 
K-matrices of the XXZ chain and give the general solution for the XXX case. 
The R-matrix of the six vertex model is given by: 



R{9) 



f 


1 











\ 







sin 7 


sinO 









sin (6+7) 


sin (6I+7) 









sin 9 


sin 7 









sin (6+7) 


sin (6I+7) 




V 











1 


/ 



(43) 
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and the general solution to the factorization equations in this model is given by 0: 

imP,X,,,()=i^""^^^'^ "^'"'V) . (44) 

^ '^'^^ \^ Asin2^ (3sm{^-e) J ^ ' 

where /3, ^, /i and A are arbitrary parameters. The associated hamiltonians to this 
K-matrix follow by the procedure used above. Defining 

K^{e)=K{9,P±,X±,fi±,^±) , (45) 
the following hamiltonians are obtained: 

N-l 

H=Y. {«+! + «+i + cosh 7 

n=l 

+ sin 7 {h-(yl - fe+cr^ + c„(Tf - c+a^ + d_a^ - c?+(t^) , (46) 
where the parameters 6±, c± and d± follow from A±, /i±, ^± and [3± as shown: 



h± = cot ^± 
2A± 



/3± sin 



± 



2/i± 



= -^^^ . (47) 
P± sm 4± 

Here /?±,^± ^ so as to have det[K-{0)] ^ and tr[fs:+(-r/)] ^ 0. 
Equation (|i6|) gives the most general choice of boundary terms compatible with integra- 
bility for the XXZ chain besides periodic and twisted b.c. In the present case one can 
twist the boundary conditions as: 



a^+i = Maf , (4J 



where M = a'' or M = e*'^'^^ <uj <2n. 
Looking to equations ([T0|) and (|l^) is now possible to see why in the elliptic case we lose a 
continuous parameter that appears in the trigonometric limit. In the trigonometric case 
d^ = d^ = and we have only the constraint of equation (|12]) which gives a continuous 



family of solutions. The same happens with eqs. (^ and (|TT|) losing again a continuous 
parameter from the elliptic case. 
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It is also interesting to see what hamiltonians are obtained from the trigonometric hmit 
of the K-matrices obtained in the preceding section. When = in ( pSf ) one obtains: 



KTAie) 



(49) 



1 G^HA sin 9 
Xha sinh 9 1 

and this is seen to correspond to solution (|^) whith /5 = 1, ^ = ±^ and /i, A = 
±A//a/2. The correspondig hamiltonians are obtained from the substitution of these 
values of the parameters in eqs. (|^), (0) and 
For solution (36) when k = one obtains: 



Ktb{9) 



(50) 



1 eXHBCOs9 
Xhb cos 9 —1 

that coresponds to solution (^) with /3 = 1, ^ = 0,7r and /i, A = ±.Xhb/'^- As 
discussed in the previous section this limit leads to a hamiltonian which includes a non- 
local coupling with the boundaries. 

It is interesting to note at this point that the trigonometric limit of the K-matrices for 
the 8- vertex model does not lead to an SUq{2) invariant hamiltonian. This is not the 
case for the free fermion 8-vertex model where the CHg{2) symmetry given by the elliptic 
K-matrices "contracts" to a Uq{gl{l,l)) symmetry in the trigonometric limit 
Let us now look for the general solution of the factorization equations in the rational limit 
of the R-matrix (S3) given by: 



m 



1 1 








1 







(e+i) 
e 1 



(e+i) 




(e+i) 




o\ 




1/ 



(51) 



The equations for the K- matrix when R is rational just follow by substituting the sine 
functions by their arguments (that is, sin(ti;) by uj) in all the equations. The number of 
independent equations is the same in the rational and trigonometric cases. (This number 
decreases going from the elliptic to the trigonometric case). Thus, the general solution is: 



ir(^,/3,A,/i,0 



m + d) f^9 

X9 /3(e - 9) 



(52) 



and using equation (|37|) one obtains: 
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N-1 



n=l 



+h^al - + c_o-^ - c+cTjv + d_a^ - d^a^ , (53) 

where we have scaled by a factor of 27 and omitted a term proportional to the identity 
operator. In this case the parameters 6±, c± and d± follow from A±, ^± and /3± as: 



^± = ^ , (54) 
P±?± 

where /3±, ^± ^ to have det[K-{())\ ^ and tr[K+{-ri)] ^ 0. 
This equation again provides the most general choice of boundary terms compatible with 
integrability for the XXX chain besides periodic and twisted b.c. 



4 Conclusions 

We have presented the general solution to the surface factorization equations for the XYZ, 
XXZ and XXX models providing in this way the most general boundary terms compatible 
with integrability. One can expect that if any kind of quantum group invariance is possible 
in the XYZ chain the necessary boundary terms will be provided by those of hamiltonian 



(41). For the XYZ chain a geneneralization of the construction for the eigenvalues and 
eigenvectors of the periodic chain remains to be done. As the hamiltonians obtained for 
the XXZ and XXX models do not commute with Jz a generalization of the Functional 
Bethe ansatz proposed by Sklyanin for open boundary conditions should be useful to 
find the eigenvalues. 

In the context of two dimensional integrable quantum field theories with boundaries it 
is interesting to solve the boundary bootstrap and cross-unitarity equations for these 
solutions. 
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